Enhanced quantum tunnelling induced by disorder 

Jean HeinrichJ*] 

Departement de Physique, B5a, Universite de Liege, 
Sart Tilman, B-J^OOO Liege, Belgium 
(Dated: May 9, 2008) 

Abstract 

We reconsider the problem of the enhancement of tunnelling of a quantum particle induced by 
disorder of a one-dimensional tunnel barrier of length L, using two different approximate analytic 
solutions of the invariant imbedding equations of wave propagation for weak disorder. The two 
solutions are complementary for the detailed understanding of important aspects of numerical 
results on disorder-enhanced tunnelling obtained recently by Kim et al. (Phys. rev. B77, 024203 
(2008)). In particular, we derive analytically the scaled wavenumber (fcL)-threshold where disorder- 
enhanced tunnelling of an incident electron first occurs, as well as the rate of variation of the 
transmittance in the limit of vanishing disorder. Both quantities are in good agreement with the 
numerical results of Kim et al. Our non-perturbative solution of the invariant imbedding equations 
allows us to show that the disorder enhances both the mean conductance and the mean resistance 
of the barrier. 
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Growing attention has been devoted in recent years to an intriguing phenomenon oc- 
curring when a quantum-mechanical particle of energy E impinges on a disordered one- 
dimensional tunnel barrier of mean height V > E. Indeed it was found that the transmis- 
sion coefficient of the barrier is increased by the effect of a weak disorder [l, [2, [J. This is 
surprising since one would expect that the additional scattering of an incoming electron by 
the disorder (e.g. random impurity potentials) would reduce the transmission rather than 
enhancing it. 

The interest in the study of the transmission of quantum particles through disordered 
one-dimensional barriers [l, 2, 2I was triggered by an earlier study of transmission of scalar 
waves through disordered one-dimensional dielectric media 4[]. In this case, the increase of 
the transmission coefficient with increasing disorder was first observed for frequencies of the 
incident wave lying in the gap of the band structure of the periodic medium in the absence 
of disorder^]. The similarity between the quantum barrier- and the dielectric medium 
problems comes from the fact that wave frequencies in the gap of the dielectric medium 
correspond to evanescent wave solutions of the wave equation, which are analogous to the 
evanescent states solutions of the Schrodinger equation for energies E < V in the potential 
barrier problem. 

Prior to the studies of transmission by disordered potential barriers in Refs. l|, [2], H], the 
author developed an extensive analytic study of resistance- and reflectance distributions in 



one-dimensional disordered tunnel barriers in the context of electron localization 



These 



were motivated by the development 



studies, based on the invariant imbedding method (3, 
of a fully probabilistic scaling theory of localization [9| in these systems where localization 
effects coexist with the familiar barrier penetration (tunnelling) effects j^l, 

The purpose of this Brief Report is first to apply the approximate solution for the reflec- 
tion amplitude of a disordered barrier obtained in Ref. 5j for deriving analytic expressions for 
mean disorder- enhanced transmission coefficients and to compare the results with previous 
works . In particular, we are interested in understanding an important new aspect 

of the transmission of a disordered barrier of length L revealed by the numerical results 
of Kihong Kim et al. jj, namely the existence of a threshold value of the scaled particle 
wavenumber kL above which transmission is enhanced by the disorder while being reduced 
by the disorder below it. 

On the other hand, it is found that the thresholds for enhanced transmission obtained 
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numerically [3J] vary significantly with the relative barrier hight V/E. In order to understand 
this behaviour, it is necessary to develop a solution of the invariant imbedding equations 
for arbitrary barriers V/E > 1 , in parallel to the solution in 5| which is specific to the case 
V/E = 2. Such a solution based on perturbation theory for weak disorder is discussed below 
and applied to study the disorder-enhanced tunnelling. The relative merits of the two types 
of solutions for describing the disorder-enhanced tunnelling are discussed towards the end 
of the paper. 

On the other hand, we close by pointing out that the non-perturbative invariant imbed- 
ding solution for a tunnelling barrier in [5j, which does lead to the phenomenon of disorder- 
enhanced tunnelling, also leads to the exponential growth of resistance and conductance of 
a disordered barrier discussed earlier by Freilikher et al.[l|. 

We consider an electron of energy E = h 2 k 2 /2m (with units such that % — m — 1) which 
impinges from the right on a random one-dimensional tunnel barrier 

V(x) = V + v(x) , (1) 

confined to the region < x < L. V denotes the mean of V(x) and v(x) is a weak gaussian 
white-noise, 

(v(x)v(x')) =£5(x-x'),(v(x)) = . (2) 
Outside the barrier the particle is described by the wavefunction 



^j( x ) = e -^-L) + r^y^-L)^ X>L ; ( 3a ) 
ip(x) =t(L)e- ikx ,x < , (3b) 

where the complex reflection and transmission coefficient amplitudes r(L) = r and t(L) = t 
are determined by the invariant imbedding equations [3] 

lk ^l = -2k 2 r(L) + V(L)(l+r(L)) 2 , (4) 
ik ^l = -k 2 t(L) + V(L)(l+r(L))t(L) . (5) 
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In I we discussed a useful approximate solution of (jl]) valid within some energy interval 
around the value E — V/2 such that for the most typical values of ([!]) one has 



2E - V ~ v(L) . (6) 

In this case the r.h.s. of (J4]) is approximately V(L)(l + r(L) 2 ) and the solution of (BJ subject 
to r(0) = is[5| 



.1 <- L 



r(L) = -itanh[- / dL'V{L')} . (7) 
k Jo 

By inserting ([71) in ([5]) we obtain the corresponding exact solution for the amplitude trans- 
mission coefficient, 

i(kL-f L dL'V{L')) 

t(L) = — - T , (8) 

cosh[I/ L ciL^(LO] 

which is valid for energies of the incident electron close to V/2. Note that ([THE]) verify 
probability conservation, 

|r(L)| 2 + |t(L)| 2 = l , (9) 

as required. 

The solution (CO) of the invariant imbedding equation (jl]) for E ~ V/2 has been further 



discussed by Haley and Erdos 



lOj in the context of the resistance, p(L), of the disordered 



barrier defined by the Landauer formula 

" (l) = T^k6f ' (10) 

These authors also developed a more general result for the Landauer resistance of a disor- 
dered barrier valid for any value of V / E. 

We first address the question of the disorder enhanced tunnelling across the potential 
barrier ([1]), for weak disorder. More precisely, by expanding |t(-£/)| 2 in ([8]) to second order in 
k~ l dL'v(L') assumed to be small and averaging over the disorder using (T5]), we get 



It follows that the sign of the effect of the disorder on the transmittance of the barrier is 
given by the sign of the factor 3 tanh 2 (VL/fc) — 1. Thus, we find that for parameters such 
that exp( ~ 2 ^ L ) < 2 — v3, the transmission coefficient of the disordered barrier is enhanced 
by the effect of weak disorder. The threshold value 

Y-k = t/vL = kL = -~ln(2- V3) ~ 0,659 , (12) 

above which the mean transmission coefficient (fTTI) increases with increasing disorder is in 
reasonable agreement with the critical value kL ~ 0, 58 at which the effect of the disorder in 
the transmittance changes sign in the results of fig. 2 of Ref . {3J . Below this threshold value 
of kL the disorder reduces the transmittance while above it the disorder enhances it. 
On the other hand, following Kim et al., we define the effective disorder parameter 

9 = ¥ ' (13) 

in terms of which we obtain from ffTTl) (with V = 2E = k 2 ) 

dg cosh 2 kL K ' V ; 

This defines the initial rate (i.e. near g = 0) of variation of the transmittance as a universal 
function of the scaled wavenumber kL, both above and below the threshold. The expression 
( TT4j) may be compared with the slopes at the origin of the transmittances as a function of 
g in fig. 2 of [2I, for various values of kL. This comparison is shown in table 1, indicating a 
quite reasonable agreement between the two sets of results. 

The enhancement of the transmittance of the tunnel barrier for weak disorder suggests, of 
course, a non-monotonic variation at larger disorder since beyond sufficiently large disorder 
the transmittance necessarily decreases to zero. This follows e.g. from the form of the 
typical transmission coefficient obtained from ([8]): 



\K L ) typical = ( ( COsh2 



2 



dL'V{L') 



2iL 2VL 
1 + e'fc 2 " cosh — - — 



n -1 



(15) 



which vanishes for £ — >• 00. The result (1151) is obtained by using the well-known formula 11 1 
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(exp 



±a / v{L')dL' 



) = exp 



a 2 £L 



(16) 



for averages over Gaussian correlated variables defined by (J2]). The non-monotonic be- 
haviour of the transmittance as a function of disorder is revealed in detail by the numerical 
calculations in Refs 0, 4]. 

We now turn to the discussion of a new simple approximate solution of the invariant 
imbedding equations, which in contrast to (JZHEJ), will be valid for arbitrary E < V. We first 
rewrite (JM5]) in terms of new amplitudes 



q(L) = e- 2lkL r(L) } s(L) = e~ lkL t{L) 



(17) 



which leads to 



ik^- = V(L) (e~ tkL + e tkL q(L)Y 



(18) 



= V(L) (1 + e™ L q(L)) s(L) 



For kL << 1 these equations reduce approximately to 



(19) 



ik^^V(L)(l + q(L)f 



(20) 



ik ^Ml ~ V (L) (1 + q{L)) s(L), kL « 1 



(21) 



and may be readily solved with the boundary conditions q(0) = and s(0) = 1. We find 



g (L) = I^y dL'V(L' 



1-4 / ^V(L' 
ik Jo 



(22) 



S (L) 



1 - 



1 



-, -1 



dL'V{L' 



, kL « 1 



(23) 



Restricting our attention to the transmittance, |t(£)| 2 = |s(L)| 2 we obtain from (0Q), ( jT71) 
and (1231) 



1 



£ (VL + j L dL>V{L>] 



(24) 
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By expanding (1241) to second order in the weak disorder v(L') in ([I]) and averaging the 
resulting expression using ([2]) we finally obtain 



(\t(L)\ 2 ) 



Q{L) 
with 



i + J,fS-i 



k 2 Q(L) \k 2 Q(L) 



(25) 



Q(L) = 1 + ^—j . (26) 

The expressions ( |25ll26|) show the existence of a wavenumber domain defined by 3 (-^r) 2 > 
1 i.e. 

kL> 7=£ • < 27 > 

where the transmittance is enhanced by the disorder. The critical scaled wavenumber thresh- 
old, -^y, may be readily compared with the critical thresholds for disorder enhanced trans- 
mission obtained numerically by Kim et al. [31] and displayed in their figures 3,4,5, successively 
for V/E = 1,5, V/E = 2 and V/E = 3. From these figures we obtain the critical values 
kL ~ 0.8 (V/E = 1,5), kL ~ 0.58 (V/E = 2) and kL ~ 0.4 (V/E = 3), whose comparison 
with the results kL = 0, 77 (V/E = 1.5), kL = 0, 577 (V/E = 2) and kL = 0.385 (V/E = 3) 
obtained from ( |271) shows remarkable agreement. 

Finally we discuss the initial slopes, d(\t(L) \ 2 )/dg\ g=0 , of the mean transmittance (1251) as 
a function of the disorder parameter ffl3l) . for various scaled incident wavenumbers kL. For 
the case VIE = 2 where numerical results for the transmittance are available in fig. 2 of 
Kim et al. [3J| , we get 

d(\t(L)\ 2 ) kL(3(kL)* - 1) 

—dj~- (l + (kLff > kL<<1 ■ (28) 
As shown in table I, the agreement between the numerical results obtained from (|28|) with 
those inferred from the results of fig. 2 of Kim et al.[3j for various kL is again quite 
reasonable. 

In conclusion, we have presented two complementary mathematical treatments demon- 
strating the existence of scaled wavenumber thresholds for the appearance of disorder- 
enhanced tunnelling of an electron and allowing furthermore to calculate analytically the 
initial rate of variation of transmittance with the disorder, as a function of the incident 
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wavenumber. Our results are in good agreement with extensive numerical calculations of 
Kim et al.^|. 

We close with a brief comparison of our results for the special energy (jS]) for large L with 
the results of Freilikher et al. l) indicating that a weak disorder increases both the mean 
resistance and the mean conductance of the tunnel barrier. The resistance (p) of the barrier 
defined by the four-probe Landauer formula is 



HL)l 



sinh 



k 



dL'V(L') 



(29) 



using ([7}l9]). For the conductance, g, the following two-probe formula 



9=W)\' , 



(30) 



is to be preferred, following well-known arguments 12j . 

From P), (USD and ([29]) we obtain at long lengths such that VL/k » 1 (with T (L) 
4exp(-2VL/k)) 



(P) 



1 2{L 



T (L) 



(31) 



and from © and (1301) 



07> = (|*(L)n~To(L)eT 



2|£ 



(32) 
1], using a 



These expressions are similar to the results obtained earlier by Freilikher et al. 
different method. 

As a final remark, we highlight the advantage of the approximate analytic solution ([THE]) 
of the invariant imbedding equations at the energy of half the tunnelling barrier height, in 
the context of the present work. On the one hand, this solution has enabled us to study the 
disorder enhanced transmittance at wavenumbers larger than those which are accessible by 
means of the perturbation analysis discussed above (see in particular table 1), which requires 
kL << 1. On the other hand, its non-perturbative character is crucial for demonstrating 
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the exponential growth of both the resistance and the conductance shown previously in [1 
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kL 


Kim et al.[3] 


Eq. dD 


Eq. d2HD 


0,3 


-0,194 


-0,205 


-0,169 


0,577 (0,58) 





-0,076 





0,6 


0,025 


-0,057 


0,019 


0,659 







0,068 


1 


0,414 


0,311 


0,25 


1,5 


0,557 


0,395 




3 


0,206 


0,058 





TABLE I: Initial {g = 0) slopes of transmittances as a function of disorder parameter g, for various 
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